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Abstract 



We report results of a systematic study of one-dimensional four-wave moving solitons 

> ! 

• i-H , in a recently proposed model of the Bragg cross-grating in planar optical waveguides with 

5_i ■ the Kerr nonlinearity; the same model applies to a fiber Bragg grating (BG) carrying 

two polarizations of light. We concentrate on the case when the system's spectrum 
contains no true bandgap, but only semi-gaps (which are gaps only with respect to one 
branch of the dispersion relation), that nevertheless support soliton families. Solely 
zero-velocity solitons were previously studied in this system, while current experiments 
cannot generate solitons with the velocity smaller than half the maximum group velocity. 
We find the semi-gaps for the moving solitons in an analytical form, and demonstrated 
that they are completely filled with (numerically found) solitons. Stability of the moving 
solitons is identified in direct simulations. The stability region strongly depends on the 
frustration parameter, which controls the difference of the present system from the usual 
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model for the single BG. A completely new situation is possible, when the velocity interval 
for stable solitons is limited not only from above, but also from below. Collisions between 
stable solitons may be both elastic and strongly inelastic. Close to their instability 
border, the solitons collide elastically only if their velocities c\ and C2 are small; however, 
collisions between more robust solitons are elastic in a strip around c\ = —02- 



1 Introduction 

Bragg gratings (BGs) are broadly used as a basis for spectral niters are other elements of 
optical telecommunication networks pp. Simultaneously, BGs offer a unique medium for the 
theoretical [2] [E] and experimental (T5] study of gap solitons (GSs), supported by the 
balance between the BG-induced dispersion and Kerr (cubic) nonlinearity dispersion . The 
latter topic was surveyed in an earlier review ^H] and recent one Matter- wave pulses 

similar to the optical BG solitons were very recently observed in Bose-Einstein condensates 

HE]. 

The BG solitons may exist not only as temporal pulses in fiber gratings, but also as spatial 
solitons in 2D (two-dimensional) waveguides ^H]- In the latter case, the BG is represented 
by a set of parallel grooves on the surface of the waveguide. A generalization, that makes it 
possible to increase the number of waves coupled by the BG, is based on a superposition of 
more than one BGs in the spatial domain. In particular, a system of three gratings which 
form a triangular lattice was considered in Ref. |TJ]. It couples three waves, giving rise to 
new types of stable GSs - for instance, "peakons". 

In recent papers [20J and [21J, a model of a cross-grating was introduced, as a superposition 
of two mutually transverse BGs written on a 2D waveguide with the Kerr nonlinearity. In 
this case, one is dealing with a system of four waves. The evolution of their amplitudes Vj, 



j = 1, 2, 3, 4, obeys the system of normal 
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Here, t is time, x and z are the coordinates in the waveguide, the Bragg reflectivity (strength) of 
one grating is normalized to be 1, and /x is a relative strength of the second grating. Coefficients 
in front of the nonlinear terms in Eqs. (0), which represent the self-phase-modulation, cross- 
phase-modulation and four-wave-mixing, correspond to the normal Kerr effect in isotropic 
materials. 

The same model finds a different interpretation in terms of a 3D photonic crystal with a 
transverse chess-board structure [2111 121] • In that case, Eqs. (JTJ), with t replaced by the third 
coordinate y, along which the medium is homogeneous, govern the spatial evolution of the 
time-independent fields. In the latter context, a model which is tantamount to the particular 
symmetric case of the system with « = 1, was earlier introduced in Ref. [22] (i n fact, this 
case is a degenerate one, as it corresponds to the gap in the model's linear spectrum collapsing 
to a point [2111 EH])- Ref. [22] was dealing with 2D solitons, treating their stability by means of 
a method which is tantamount to the known Vakhitov-Kolokolov (VK) criterion [23] • In fact, 
this method is not sufficient for the verification of the stability, and our direct simulations (to 
be reported elsewhere) demonstrate that all the 2D solitons, both in the general system (JT|) 
and in its special case with /j = 1 that was considered in Ref. [22], are unstable. 

A simpler possibility is to consider ID solitons, which are generated by the substitution 
Vj(t,x,z) = e tkz Uj(t,x) in Eqs. (P), 

,du\ k i du\ 

l ^t~V2 Ul+ V2~dx~ 
du 2 k i du 2 

,du 3 k i du 3 

2—^ ^u-i — — — 

dt y/2 V2 dx 

,du4 k i <9«4 

Note that Eqs. ([2]) conserve the Hamiltonian, momentum, and the norm (which is frequently 
called energy, in the applications to optics), 

4 P + OO 

E = J2 luj^fdx. (3) 

3=1 J -°° 

The model based on Eqs. ([2]) is of special interest, as previously considered multi-wave 
GS models, such as the three- wave one in Ref. [21], did not contain a parameter similar to k, 
which accounts for frustration, i.e., symmetry breaking in the wave pairs (^1,^2) and (^3,^4) 
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in Eqs. (J2J). In fact, the frustration makes the present four-wave distinct from the usual 
two- wave BG model. 

The same system (J2J) also finds an interpretation in terms of the BG in optical fibers. 
In that case, the fields and 1*3,4 in Eqs. (J2J) represent pairs of right- and left-travelling 
waves with two orthogonal circular polarizations, \x is the Bragg-reflectivity strength, while 
the linear-coupling terms with the coefficient 1 account for linear mixing between the circular 
polarizations due to an elliptic deformation of the fiber's core, and the terms proportional to 
k take into regard birefringence due to a twist of the fiber. 

Looking for solutions to the linearized version of Eqs. (J2J) as Uj(x,t) ~ exp (iXx — iut), 
one arrives at a dispersion relation, 

W 2 = i(V2/i 2 + A 2 ±v / 2 + P) 2 , (4) 

which generates a bandgap in the spectrum, |u;| < uj e dgc = fi — yl + (P/2), provided that 
/i 2 > 1 + (k 2 /2). Besides the gap, the two branches of the dispersion relation (p£J) give rise to 
semi-gaps, max j/x — y/l + (k 2 /2), j < |u;| < \i + y/l + (k 2 /2), which are also found in the 
case of 

< 1 + (fc 2 /2) , (5) 

when the full gap is absent. It was recently shown, in a general form, that families of BG 
solitons may exist in semi-gaps even if the full gap is absent [21]. 

Quiescent (zero- velocity) solitons in Eqs. (J2J), of the form Uj(x,t) = Uj (x) exp (— kot) , 
were investigated in detail in Ref. |20J [solitons generated by Eqs. (JIJ in the spatial domain, 
corresponding to Uj(t, z, x) = exp (—iut) Uj(z, x) with fixed lu, were studied in Ref. jZI]]- They 
obey a reduction U3 = —U*, U4 = —U%, which is not valid for moving solitons, see below. 
If k = 0, the solitons fall into three categories: symmetric and anti-symmetric ones, with 
Ui = ±[^2, and more general asymmetric solitons. The latter ones, and the solitons obtained 
by their continuation to k 7^ 0, are always subject to a weak instability. The symmetric 
and anti-symmetric solitons, as well as those stemming from them at k 7^ 0, were found, 
respectively, to be stable chiefly (but not exclusively) in the full gap (if any), and in the 
semi-gaps outside the full gap. 

Moving solitons are to be looked for solitons in the form 

Uj(x, t) = exp {—kot) Uj(£), £ = x — ct. (6) 
4 



Straightforward consideration of the underlying system (J2J) shows that moving solitons may 
exist for |c| < c max = l/v2 • On the other hand, quiescent BG solitons have never 
been observed in the experiment, the minimum velocity at which they could be created being, 
roughly, c max /2 [15J (a possibility to diminish the velocity through collisions between the 
moving solitons was proposed in Ref. [25] )• In Ref. [2H], moving solitons in the present model 
were considered, in a brief form, inside the full gap. It was found that they loose their stability 
at quite small values of the velocity, ~ 0.2c max . Our aim is to study the stability of the moving 
solitons in a systematic way, as well as collisions between them. We will consider moving 
solitons in the most interesting case, when the full bandgap does not exist, but, nevertheless, 
solitons may be found in the semi-gaps (see above). To the best of our knowledge, the moving 
solitons in this case has never been considered before. 

2 Existence and stability of the moving solitons 

To look for solitons in the form (JHJ), one can rewrite the underlying equations (J2J) in terms of 
the variables (t, £). The linearization of the equations gives rise to a dispersion relation, which 
differs from (0J by the replacement uj — > uj + cA. In the case when the spectrum contains 
no full gap, straightforward analysis yields an equation to for the values of uj at edges of the 
semi-gaps, where we aim to find the solitons: 



It immediately follows from Eq. (JJJ) that no semi-gaps, hence no solitons either, are possible 



Stationary equations, which are obtained by the substitution of Eq. (JBJ) into Eqs. (J2J), were 
solved by means of the shooting method. The numerical solution shows that the semi-gaps 
are, indeed, completely filled by solitons. Note that the moving solitons cannot be generated 
automatically from quiescent ones, as Eqs. (J2J) bear no Galilean or Lorentzian invariance. 

Stability of the solitons was identified in direct simulations of Eqs. (J2J) for perturbed 
solutions, by means of the split-step Fourier-transform and finite-difference techniques, both 
methods yielding identical results. Typical examples of the stable moving solitons, found 
inside the semi-gaps, are displayed in Fig. [2] 




(7) 



for |c| > l/y/2. Typical examples of the semi-gap in the plane (c,uj) are shown in Fig. ^ 
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Fig. 1: The solid curves are edges of the semi-gap, determined in the analytical form by Eq. © with 
the lower sign in the left-hand side [a counterpart of this semi-gap, corresponding to the upper sign in 
Eq. 171), is located at lo < 0]. The region between the solid curves (the semi-gap proper) is completely 
filled by the solitons, the dotted curve being the border between stable and unstable ones, as identified in 
direct simulations. Shown are typical examples for the cases of zero frustration, k = 0,/i = 0.5 (a), small 
frustration, k = 0.005,/i = 1 (b), and large frustration, k = = 1 (c). Three double arrow lines in (b) 
are cross sections of the stability region for u = const, which correspond to Eqs. 
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Fig. 2: Generic examples of stable moving solitons in the form of Eq. found by the shooting method. 
Shown are real and imaginary parts of the four fields. The panel (a) and (b) represent, respectively, the 
cases of zero (k = 0, [i = 0.5, c = 0.5, u = 1) and finite (k = 1, = 1, c = 0.2, uj = 2.13) frustration. 
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As well as in the case of zero- velocity solitons (201, the frustration parameter k strongly 
affects the solitons' shape and stability: if k = 0, the stability area is largest, and the stable 
solitons obey the reduction U\$ = —U%^. As is seen in Figs. [2] and [H nonzero frustration 
breaks the latter relations, and leads to strong shrinkage of the stability region. 

A notable feature is that, in the case of small but finite k, the stability interval for a fixed 
value of uj may be limited not only by a maximum value of the soliton's velocity, but also by a 
minimum one. For example, three cross sections shown by arrows in Fig. Wlb) (for k = 0.005 
and \i = 1) correspond to the stability intervals 

< |c| < 0.3082 for u = 1.9; 
0.28 < |c| < 0.5049 for to = 1.7, 0.6595 < |c| < 0.70356 for to = 1.1, (8) 

the two latter ones being limited by a minimum value of c. The fact that the stability of the 
solitons may require a finite velocity is a completely new one. Nothing similar occurs in the 
model supporting the usual two-wave GSs [I]. 

The stability of solitons is not sensitive to the value of the relative Bragg reflectivity 
unlike the frustration k. In particular, the variation of fi does not essentially change the shape 
of the stability region displayed in Fig. ^ nor the shape of the solitons shown in Fig. [21 

Note that stability of the solitons against non-oscillatory perturbations (those correspond- 
ing to real stability eigenvalues) is provided, at fixed c, by the above-mentioned VK criterion, 
dE/du < [2B] [recall the soliton's energy E is defined as per Eq. (J3J)]. The applicability 
of the VK criterion to the present model was not rigorously proven, but in Ref. [221 ^ was 
actually derived for the particular case of fi — 1. Inspection of numerical data shows that all 
the solitons satisfy this criterion; however, it does not preclude instability against oscillatory 
perturbations (those corresponding to pairs or quartets of complex instability eigenvalues), 
that is why a part of the solitons are unstable. Further simulations demonstrate that the de- 
velopment of the instability always leads to complete destruction of the soliton, see a typical 
example in Fig. El 

3 Collisions between moving solitons 

The next step in the study of the moving solitons is to consider collisions between them, 
which is interesting in its own right, and for various applications too [22]. As Eqs. © are 

8 



Fig. 3: A typical example of the destruction of an unstable soliton, observed in direct simulations of Eqs. 

Here (as well as in Figs. 4 and 5 below) shown is the evolution of the field \ui\; in other components, 
the evolution is quite similar. The parameters are k = 0, \i = 0.5, c = 0.4, uj = 0.6. The figure uses the 
original coordinate x, rather than £ from Eq. 

not integrable, an issue is whether collisions between stable solitons will be elastic or not. We 
stress that neither the stability of the solitons of the semi-gap type, considered in this paper, 
nor collisions between them were considered in earlier works. 

Exploring the stability region of the moving solitons (see Fig. [TJ, we have found that 
collisions between them may be both completely elastic, leading to no other effect but finite 
shifts of the colliding solitons, and strongly inelastic, resulting in a partial destruction of the 
solitons, and generation of large amounts of radiation. Typical examples of these outcomes 
are shown in Figs. 0] and EJ The interaction between the elastically colliding solitons has the 
character of the bounce from each other. 

The outcome of the collision strongly depends on the solitons' velocities ci^. Therefore, 
the results can be naturally summarized in the form of a region in the plane {01,0%) where 
the collisions are elastic; at the border of the region, the transition between the two types of 
the collision is found to be quite sharp. In turn, the size of the region of the nondestructive 
(elastic) collisions depends on proximity of the individual solitons to the instability border in 
their parameter planes, see Fig. ^ 

First, in Fig. El we present these results for the zero-frustration case (k = 0), when the 
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Fig. 4: A typical example of the elastic collision between two identical stable solitons moving with the 
velocities C12 = ±0.3, the other parameters being k = = l,u = 2.08. 




Fig. 5: A typical example of the inelastic collision between two stable solitons moving with the velocities 
ci = —0.1 and C2 = +0.25. The other parameters are k = 0,fi = 0.5, uo = 1. 
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stability region in Fig. HJa) is quite large. As it follows from the inspection of Fig. QJa), 
the difference between the three cases displayed in Fig. |B] is that in the case of u = 1 (a) 
the colliding solitons are close to the instability border, while in the other cases, u — 1.2 
(b) and to = 1.4 (c), they are located deeper in the stability region (in the latter case, the 
soliton is actually close to the edge of the semi-gap, i.e., the existence border of the solitons). 
Accordingly, there is a notable difference in the shape of the elasticity region in the {c\,C2) 
plane: in the first case, the collisions may be elastic only if the velocities are small, while in the 
other cases, the solitons with large but nearly opposite velocities, c\ ~ — c 2 , collide elastically 
even if the velocities are large. 

As it was stressed above, in the case of finite frustration, the velocity interval corresponding 
to the stable solitons may be limited not only from above, but also from below. The collisions 
between unstable solitons being definitely destructive, the latter fact strongly affects the region 
of elastic (nondestructive) collisions in the (ci,c 2 ) plane for k ^ 0, as is seen in Fig. (J7J). 

As it was mentioned above, in the currently feasible experiments with BG solitons, their 
velocity cannot be made smaller than, approximately, half the maximum group velocity, which 
is 1/ \/2 in the present notation. For this reason, it is quite important that the elasticity regions 
shown in Figs. EIb,c) and[7{a,b) extend to large values of c. 

Lastly, it is necessary to mention that, as well as the stability of individual solitons, the 
outcome of the collisions is not sensitive to the relative BG reflectivity fi, on the contrary to 
the strong dependence on the frustration k. In fact, the elasticity regions show almost no 
dependence on /i, therefore the examples shown in Figs. |H]and[7|are quite generic ones. 

4 Conclusion 

In this paper, we have presented results of a systematic study of the stability and collisions 
of moving four-wave solitons in the recently introduced model of the Bragg cross-grating in 
planar optical waveguides with the Kerr nonlinearity. We have focused on the case when the 
system's linear spectrum does not contain any true bandgap, but rather features semi-gaps, 
in which soliton families may be found. In recent works, only zero-velocity solitons in this 
four-wave system have been studied. That case is far from the current experiments, which 
generate solitons whose velocity exceed, roughly, half the maximum group velocity. We have 
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Fig. 6: Areas of elastic and inelastic collisions between solitons in the plane of their velocities, (ei,C2) 
(the values up to the limiting value, c max = l/v2 « 0.71, are shown). The fixed parameters are k = 
and \i = 0.5. The frequency is u = 1.0 (a), 1.2 (b), 1.4 (c). The horizontal and vertical dashed lines 
limit the velocity intervals of the stable solitons, as per Fig. U 
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Fig. 7: The same as in Fig. © for the case of nonzero frustration, k = 0.005 and /x = 1. The fixed 
frequency is = 1.1 (a), 1.7 (b), 1.9 (c). 
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found the semi-gaps for the moving solitons in the analytical form, and demonstrated that 
numerically found solitons completely fill the semi-gaps. Stability of the moving solitons was 
determined by direct simulations, revealing a stability border which strongly depends on the 
frustration parameter. It is interesting that a previously unknown situation is possible in the 
present model, when the velocity interval for stable solitons is limited not only from above, 
but also from below. 

Collisions between stable solitons may be both elastic, so that they bounce from each 
other, and strongly inelastic, leading to a partial destruction of the solitons, with a sharp 
border between the two outcomes. A region of elastic collisions was identified, with a well- 
pronounced peculiarity: if the individual solitons are close to the stability border, they collide 
elastically only if their velocities c\ and c 2 are small; however, if the solitons are more stable, 
the collision remains elastic in a strip around the line c\ = — c 2 . As well as the stability of 
individual solitons, the elasticity region strongly depends on the frustration parameter, but 
is insensitive to the relative reflectivity of the two Bragg gratings which constitute the cross 
configuration. 
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